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In the coordinate representation of thermofield dynamics, we investigate the thermalized dis- 
placed squeezed thermal state which involves two temperatures successively. We give the wavefunc- 
tion and the matrix element of the density operator at any time, and accordingly calculate some 
quantities related to the position, momentum and particle number operator, special cases of which 
are consistent with the results in the literature. The two temperatures have different correlations 
with the squeeze and coherence components. Moreover, different from the properties of the posi- 
tion and momentum, the average value and variance of the particle number operator as well as the 
second-order correlation function are time-independent. 
Introducing finite temperature effects into squeezed states is important, because a squeezed state can possesses 
minimum uncertainty and squeezability and accordingly technological applicability |J and, on the other hand, a 
finite-temperature influence on it is inevitable. This problem has received extensive investigations [fq-|llfl. Analyzing 
these investigations, Ref. H,|| divided the squeezed states with finite temperature effects into the thermalized squeezed 
states and the squeezed thermal states • The thermalized squeezed states and the squeezed thermal states 

£s^j | are physically distinct states, albeit they are turned into each other by some parameter transformation ||. Each 
of these two states has its several possible representations, and Ref. Q gave a detailed discussion about them and 
— i ' elucidated their physical interpretations. From Ref. it is not difficult to understand that the squeezed thermal states 
correspond to the output from a squeezed device whose input is a thermal chaotic state with Bose-Einstein distribution, 
while the thermalized squeezed states are prepared by thermalizing a squeezed state provided the thermalizing source 
is such a device that it can bring a vacuum state into a thermal chaotic state. Both the thermalized squeezed state and 
the squeezed thermal state are involved in one thermal source only. In this paper, we intend to consider yet another 
different state which we call the thermalized displaced squeezed thermal state (TDSTS). This state can be prepared 
ON i by thermalizing a displaced squeezed thermal state, and obviously contains the thermalized displaced squeezed state 
ON " (TDSS) and the displaced squeezed thermal state (DSTS) as its special cases. The TDSTS is involved in two thermal 
sources successively and perhaps is more practical than the TDSS and the DSTS, because, in general, there exist 
thermal noises for the input and output of a displacement-squeeze device which generates a displaced squeezed state, 
and furthermore the thermal noises have perhaps different temperatures for the input and output, respectively. 

In this paper, within the framework of thermofield dynamics [[l2||l3| l, we shall first give the definition of the ther- 
malized displaced squeezed state, then give the time-dependent wavefunction and calculate the matrix elements of the 
density operator in the coordinate representation. From the density matrix elements, the probability density, average 
value, variances and the entropy uncertainty relation of the position and momentum will be calculated. Finally, 
, we shall also give the average value and variance of the particle number operator and the second-order correlation 
J^j ' function. 

This paper is based on the one-dimensional quantum oscillator with a mass m and a constant angular frequency lo 

■ whose Hamiltonian is 

_ _ i 

H = —p 2 + -muj 2 x 2 = (a^a + -)Tiuj , (1) 
2m 2 2 

where x is the position operator, p = —ih4z = —ihd x the momentum operator in the coordinate representation, and 

a = (ip + mux) , a) — {~W + muix) (2) 

\2mTiuj \2mhio 

are the annihilation and creation operators, respectively. Nevertheless, taking the mass as unit in the formulae of this 
paper, one can get the results which are usable for a one-mode electromagnetic field with the same frequency. 

As is well known, a displaced squeezed state can be constructed by the squeeze and displacement operators succes- 
sively acting on the ground state |0 > of the oscillator Eq.(l) mM. The displacement operator 
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D(a)=e aa aa , (3) 

with a = {a,\ + ia.2) = |a|e 17 any complex number, corresponds to an ideal displacement device ( the symbol '| • • • |' 
represents the module of a complex number ), and the squeeze operator 

S(z) = exp{-i(z*aa - zaV)} , (4) 



with z = z\ + iz2 — re 1 ^ any complex number, corresponds to an ideal squeeze device. ( In Eq.(4), the minus "— " 
before "i" is sometimes replaced by the plus "+" in the literature, which gives rise to no essential differences. ) The 
action of the displacement operator following squeeze operator on the ground state will yield the displaced squeezed 
state D(a)S(z)\0 > ( we call the corresponding device the displacement-squeeze device ), and a squeezed displaced 
state S(z)D(a)\0 > can be also constructed by the action of the squeeze operator following the displacement operator. 
Note that through a parameter transformation the states D(a)S(z)\0 > and S(z)D(a)\0 > can be turned into each 
other, because one has [|l5| 

S(z)D{a) = D(a cosh(r) + a V sinh(r))5(z) . (5) 

In order to consider finite temperature effects, thermofield dynamics [[l2| introduces a copy of the physical oscillator 
Eq.(l) ( called the tilde oscillator ) 

H = ^—p 2 + \mtu 2 x 2 = (a) a + \)Tiuj , (6) 
2m 2 2 

according to the tilde "conjugation": CO = C*0 Here, C is any coefficient appeared in expressions of quantities 
for the physical system, O any operator, the superscript * means complex conjugation, and O represents the corre- 
sponding operator for the tilde system. Based on the physical and tilde oscillators, thermofield dynamics manufactures 
a thermal operator T(9) 

Q 

T{9) = cxp{-6»(/3)(aa - a)a<)} = exp{i-(xp - xp)} (7) 

with 

tanh[6>(/?)] = e~ (ihwj2 

and /3 = -^p- Here, k}, is the Boltzmann constant and T the temperature. The thermal operator is invariant under 

the tildian conjugation, i.e., T(0) = T{9). Letting T{9) act on the direct product of the physical ground state |0 > 
and the tilde ground state |0 >, one can manufacture a thermal vacuum. Notice that any physical operator commutes 
with any tilde operator, physical operators act on physical states only, and similarly tilde operators on tildian states 
only. Consequently thermal-vacuum average values in thermofield dynamics agree with canonical ensemble average 
values in statistical mechanics jl2| . When the thermal operator acts on the direct product of the free elctromagnetic 
field vacuum and its tilde counterpart, one can get a thermal chaotic state which describes a thermal chaotic light 
with a Bose-Einstein distribution at the temperature T . That is to say, generally, the thermal operator T(6) can 
correspond to the action of a thermal source. 

Now we can give the definition of the TDSTS. Since T(9) is involved in the physical and tildian operators simul- 
taneously, we shall use both the physical operators D(a), S(z) and their tildian counterparts D(a), S(z) when a 
finite temperature effect is introduced into a displaced squeezed state ( The tildian operators D{a) and S{z) were not 
adopted in Ref. Q). Because T{9) commutes with S(z)S(z) || and does not commute with D(a)D(a), we can define 
a state which involves two thermal sources with different temperatures successively. Thus, we have the following 
TDSTS 

\f32,a,z,/3 1> 0>=r(e 2 )D(a)D(a)S(z)S(z)T(ei)\0> |0 > , (8) 

where (3\, 9\ and [32, #2 correspond to those at the temperatures T\ and T2, respectively. Furthermore, the time 
evolution of the TDSTS can be considered by the evolution operator U(t) — exp{ — i(H — H)} acting on the TDSTS, 
i.e., 

\t,j3 2 ,a,z,Pi,Q >= U(t)\p 2 ,a,z, Pi,0 > (9) 
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with t the time. Obviously, the DSTS §||§ is the TDSTS for T 2 = 0, the TDSS for Ti = 0, and the thermalized 
coherent thermal state for z = 0. The TDSTS Eg. (8) with proper parameter constraints and transformation can be 
reduced into almost all cases discussed in Refs. p|-pl| except for the displaced thermalized squeezed states defined by 
a density matrix in Ref. B. The TDSTS is physically meaningful. For example, when a detector examines the output 
from the displacement-squeeze device whose input is a thermal chaotic state at the temperature Ti, a TDSTS with 
Ti and T2 effects will be recognized by the detector provided the detector's added noise is identified with a thermal 
chaotic state at the temperature T2. A thermal chaotic state can be regarded as the ground state with a thermal 
noise at some temperature. Thus, the TDSTS involves two thermal noises successively. For convenience sake, we call 
the thermal noise at T2 the detector thermal noise and the thermal noise at T\ the input thermal noise. 

In the coordinate representation, the displacement, squeeze and thermal operators all can be unentangled p6|,pT| . 
Moreover, if exploiting the thermal coordinate representation introduced in Ref. (T^j, one can easily uncntangle 
the evolution operator U(t). Hence one can obtain the explicit expession of the time-dependent wavefunction < 
x,x\t,P 2 ,a,z,Pi,0 >. Alternatively, noticing that D(a)D(a)T(6) = T(6)D(a(cosh(6) - sinh(6))))T)(a(cosh(6>) - 
sinh(6>))) |lq [, one has \t,02, a, z,/3i, >= \t, jj 2 , a(cosh(fli) — sinh(#i)), z, > and can obtain easily the time- 
dependent wavefunction with the help of the result in Ref. p7| . Taking n = 0, 9 = 61 + 9 2 and the replacement 



a — ► a(cosh(6>i) — sinh(0i)) in Eq.(44) of Ref. |17|, one can read 



<x,x\t,0 2 ,a,z,p 1 ,O> =(— ) — exp|--- — } 

• exp I — -t^-Gi ( x cosh(0) — x sinh(6)) z 



2% 

I mm 

-2 



y / ^G 2 (cosh(6<i) - sinh(6» 1 ))(xcosh(e) - x sinh(9))} 

' exp { ~ 1>K G * cosh ( e ) ~ * sinh ( )) 2 



+2W^G^(cosh(6'i)-sinh(6i 1 ))(icosh(e)-a;sinh(e))} , (10) 



where 



T\ = cosh(r) + sinh(r) cos(0) + isinh(r) sin(</>), 
^_ 1 — % sinh(2r) sin(0) 

2 cosh(2r) + sinh(2r) cos(</>) ' 
9 = 6»i + 9 2 , B = cos(wi) + 1T2 sin(wt), 

Ti cos(ujt) + i sin(wi) + ia 2 

( 1 1 — Gro — , 

S B 

and 

Q = [-^2 cos(wi)a 2 + 2^2 sin(w)aia2 + i sin(wi)a2](cosh(#i) — sinh(^i)) 2 . 

Eq.(10) is the time-dependent wavefunction of the TDSTS in the coordinate representation, and all information of 
the TDSTS can be extracted from it. 

With the help of Eq.(10), a straightforward calculation yields the density matrix element p x >. x (t) of the position 
for the TDSTS 



(t) = I < x, x\t, P2, ot, z,Pi, >< 0, Pi, z, a, 02, t\x , x > dx 

^shtm) exp ( 2 cosh(2e) coth (— - ^ ) 



r muj 1 

• eX Pi " 7i n|2, 



47i |J r i J B| 2 cosh(26) 



cosh(20) 



27i / , ./^Tiw 



x - x' - \/ ; /»^x \/coth(^— )(Q f 2 - G2) 

w muj cosh(26) w 



2 



47i l^iBI 

.^(Gx-GDC^-x' 2 )} (11) 
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with A = cos(wt) + isin(tjt). This density matrix is complex, Hermitian and time-dependent. When z — Eq.(ll) 
gives the density matrix of the position for a thcrmalized coherent thermal states T(02)D(a)D(a)T(6)\O > |0 > which 
contains the thermalized coherent state and the coherent thermal state as its special cases. At the initial time t = 
and T2 = 0, Eq.(ll) is reduced into the density matrix element of the position for the DSTS 



Imoj 1/1 r LFiT r . l2 l 

Px'.x = \/ % , 4/ , , „ s exp<^ __ — J- 2 _ J" )ai + «2a 2 f 

7rft |.Fi| y cosh(20i) L2cosh(2#i) J 

71C<J 1 

4ft _ |.7 : i| 2 cosh(26>i) 



f 1 r , / 2ft n i 2 mw , _ __,, w 9 ,o s 

I 4ft LFi 2 cosh(26»i ) L V mw J 4ft v ^ v ; 



muj cosh(26'i) 
4ft 1^1 1 2 



x — a; 



raw cosh(26> 



-((^ 2 -^ 2 *)a 1 +i2a 2 )l } , (12) 
1) J J 



which is identical with Eq.(6.5a) in Ref. §(1993). Of course, setting Ti = 0, one can obtain the density matrix 
element of the TDSS. Noticing 

cosh(2B) = coth( — - — ) coth( — - — ) + cosech( — - — )cosech( — - — ) (13) 

and cosh(6>) = coth(^), one can see that the differences both among the density matrices of the TDSTS, the 
DSTS as well as the TDSS and among the finite temperature influences on these states consist in the appearance 
or disappearance of the four factors cosh(26), cosh(2#i), cosh(20 2 ) and coth(^M^). In the expression of the density 
matrices, the factors cosh(26), cosh(26>i) and coth(^£) appear for the TDSTS, the factors cosh(20 2 ) and coth(^^) 
for the TDSS, and for the DSTS the factor cosh(26»i) replaces cosh(26). 

Taking x' = x in Eq.(ll), we have the probability density of the position for the TDSTS 



Imu! 1 / 1 

Px,x{t) 



?rft \FiB\V cosh(29) 



r mui \ 
•exp^ 



ft |J c 'iB| 2 cosh(2e) 



2muj 



coth(^ r )(- + -)] }. (14) 



This is a Gaussian distribution, and from it one has easily the average value of the position 



< x > ee I xp x , x dx = ; /_Wcoth(— )( 3 + - 



= \/— ycoth(-2— )|a|cos(a;i-7) (15) 
" muj V 4 



and the position variance 



(Ax) 2 = <x 2 >~<x> 2 = a| f li?l cosh(26) 

cosh(2r) + sinh(2r) cos(2wi - 0)] cosh(26) . (16) 



2mu> 



Here and after, "< • • • >" denotes "< 0, /3i, z, a, /3 2 , t\ ■ ■ ■ \t, 02, a, z, fti, >". 

Furthermore, exploiting Eq.(ll), one can calculate the probability density of the momentum for the TDSTS 



/OO 1 / 
1 PX PX ill 
-^-^ exp{i— i — \p x >. x dxdx 



1 \FxB\- 



irmfiLO 
■exp\ - 



,Gj-G* 

\FiB\ 2 



1 - |.FiB|4(£lzS.)2 V cosh(26) 



mftw[l - l^iBl^^^) 2 ] cosh(28) 
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co th(— — ) (- + — ) - (G 2 - G 2/ 




2 



1)21 cosh(20) 



p + 1 



coth( — )(---)] }. (17) 
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(Noting that in this paper, some symbols, for example, p, p, should be understood as either operators or the corre- 
sponding eigenvalues of the operators according to the context where they appear.) The probability density of the 
monmentum is also Gaussian. The average value and variance of momentum are 



and 



< P >= -*V — coth (^7-)(]4 -#;) = -\J 2m ^ coth(^) ainM - 7 ) (18) 



(Ap) 2 = < p 2 > - < p > 2 = [ - 2 ' J cosh(29) 

= [cosh(2r) - sinh(2r) cos(2tjt - 0)] cosh(26) , (19) 



respectively. 

When T 2 = and t = 0, Eqs.(14) and (17) are consistent with Eqs.(6.6) in Ref. [§(1993). For the case T x = and 
t = with both m and fi. unit, if taking the replacement a — ► [acosh(r) + a*e^sinh(r)] ( Eq.(5) ) in Eqs.(15),(16),(18) 
and (19), one can obtain correspondingly the self-tildian results in Ref. 0(1989). Eqs.(15) and (18) indicate that 
for the TDSTS the average values of the position and the momentum follows the harmonic motion of a classical 

oscillator. The detector thermal noise widens the amplitudes of the harmonic motion by the factor J coth(^|^), but 

the input thermal noise takes no effects on < x > and < p >. However, from the factor cosh(29) of Eqs.(16) and 
(19), the input thermal noise enlarges the variances of the position and momentum at the same way as the detector 
thermal noise. It is worth notice that the joint enlargement factor cosh(28) of the variances by the input and detector 
thermal noises is greater than the product of the corresponding enlargement factors for the TDSS and DSTS, namely, 
cosh(2Q) > coth(^) coth(^) ( see Eq.(13) ) . 

According to Eqs.(16) and (19), the product of the standard deviations of the position and momentum is 



(Ax)(Ap) = - cosh(2e)ycosh 2 (2r) - sinli (2r) cos 2 (2wi - <f>) (20) 

which is consistent with the thermal Heisenberg uncertainty principle fH| . Obviously, at any time, neither is the TD- 
STS a minimum Heisenberg uncertainty state, nor a minimum thermal Heisenberg uncertainty state. This is different 
from the TDSS and the DSTS which can be a minimum thermal Heisenberg uncertainty states when cos 2 (2o;i — <p) = 1. 

From Eqs.(14) and (17), one can also obtain the sum of the information entropies of the position and the momentum 
in the TDSTS @ 

/OO pOO 
Px,x{t) \n[p x . x (t)]dx - / p p . p (t)\a[p p . p (t)]dp 
-OO w — OO 

= ln[7re7i] + ln[cosh(29)] + - ln[cosh 2 (2r) - sinh 2 (2r) cos 2 (2wi - <j>)} (21) 

which is in conformity with the thermal information-entropy uncertainty relation fll9|| . Due to cosh(20) > cosh(2#2), 
it is impossible that the TDSTS is a minimum entropy uncertainty state. 

The squeeze effect can be discussed by the rotated quadrature phase operators (the definitions here are slightly 
different from those in both Ref. [|| and Ref. [Q(the book)) 

Y 1 = hae-^+a^e i ^), Y 2 = ^(ae"** - Je*) , (22) 
I 1% 
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which give the quadrature phase operators X\ = y^2jf x and X 2 = y 2r ^ huJ P wnen the rotated angle (f = 0. A 
straightforward calculation yields 

(AYi) 2 =< Y? > - < Yi > 2 = i cosh(26)[cosh(2r) + sinh(2r) cos(2wt + 2y? - 0)] (23) 

and 

(Ay 2 ) 2 =< K 2 > - < Y 2 > 2 = j cosh(26)[cosh(2r) - sinh(2r) cos(2wi + 2^ - 0)] . (24) 

The last two equations indicate that for the TDSTS, it is always possible to attenuate either AYi or AY 2 at any value 
of (j) and any time t. At t = 0, Eqs.(23) and (24) are identical to Eq.(6.9) in Ref. || (1993) when tp = 0, T 2 = 0, and 
to Eq.(3.7) in Ref. § when T 2 = and (0 - 2ip) = n. 

For the probability densities Eqs.(14) and (17), the moment generating functions of the position and the momentum 
are 

/OO 1 
e Xx Px , x (t)dx = exp{A < x > +-A 2 (Ax) 2 } (25) 
-oo " 



and 



f°° 1 
M P (X) = J e x P Pp Jt)dp = ex P {\<p>+-\ 2 (Ap) 2 } , (26) 



respectively. Hence one can calculate the nth moment of the position as per 

dM x {\) 



< x" >- 



and the nth moment of the momentum as per 

dM P {X) 



< p" >- 



dX 



A=0 



A=0 



(27) 



(28) 



They can be used to discuss the higher-order squeeze effects for the TDSTS j2l| Q(1991). 

We have discussed the properties related to the position and the momentum. Next, we will give the average value 
and variance of the particle number operator n = {^^jK 1 ^) 2 — ®x\ ~ h} an d the second-order correlation function. 
The average value of n at any time t is 

< n >= 1 cosh(26) cosh(2r) - i + coth(^ 2 |^)|a| 2 . (29) 

Thus, the input thermal noise is correlated only with the squeeze, while the detector thermal noise take effects both 
on the squeeze and on the coherence. When T 2 = 0, Eq.(29) is Eq.(3.9) in Ref. |](1993) and Eq.(3.3) in Ref. [§. 
Exploiting Eqs.(27) and (28), one can calculate the variance of n at any time and obtain 

(An) 2 = < n 2 > - < n > 2 = -- + - cosh 2 (26) cosh(4r) 

+ cosh(2e)coth(^ 2 |^)|a| 2 [cosh(2r) + sinh(2r) cos(0 - 2 7 )] . (30) 

Consequently one has the second-order correlation function at any time t 

9^(0)= <n2> - < ? n> =2 + 



< n > i 

\ cosh 2 (20) sinh 2 (2r) - coth 2 (^)|a| 4 + cosh(28) coth(^)|a| 2 sinh(2r) cos(0 - 2 7 ) 



\\ cosh(29) cosh(2r) - \ + coth(^)|a 



212 



(31) 



When T 2 = the last equation is identical to Eq.(3.10) in Ref. § (1993). From Eqs.(29),(30) and (31), one sees 
that the average value and variance of the particle number operator and the second-order correlation function are 
time-independent. 
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In conclusion, this paper has investigated the TDSTS which is a generalization of the TDSS and the DSTS. In 
the coordinate representation of thermofield dynamics, we give the wavefunction and the density matrix with time 
evolution, and calculate some quantities related to the position, momentum and particle number at any time. Our 
results indicates that the quantities related to the position and the momentum are time-dependent, but the average 
value and variance of the particle number operator and the second-order correlation function are time- independent. For 
the influence of the temperature on the displaced squeezed state, the TDSTS possesses the features both of the TDSS 
and of the DSTS, but the finite-temperature effects of the TDSTS are not a simple sum or product of those of the TDSS 
and the DSTS. The input thermal noise at T\ influences only the squeeze via its contribution to the factor cosh(20), 
while the detector thermal noise at T 2 has an additional effect on the coherence component by the factor coth( ^" ). 
The TDSTS contains various thermal squeezed states in the literature as its special cases, except for the density- 
matrix-defined displaced thermalized squeezed state in Ref. For the definition Eq.(8), if inserting the thermal 
operators T{9^), T(#4) and T(#s) between D(a), D(a), S(z) and S(z), respectively, without self-tildian requirement 
(as Ref. @(1989) did), one can obtain a most generalized displaced squeezed state with finite-temperature effects. 
However, we do not know if the most generalized state can be reduced the density-matrix-defined displaced thermalized 
squeezed state in Ref. Q. Finally, we want to point out that, if there are M input thermal noises and N detector 
thermal noises and the corresponding temperatures are T± t i, Ti,2, • • • , Ti,m and T^, • • • , T^jv, then the physical 
quantities of such a generalized TDSTS can be given by taking 8 = + T 1;2 + • ■ ■ + + 7*2,1 + ^2,2 + ■ • ■ + T 2 .n 
and T 2 = T2.1 + T2,2 + ■ ■ ■ + T^jv m the results of this paper. 
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